/* Create Figure 1 from Model 1 */
* Note the blindscheme add-on is needed to recreate the format of the figure *

/* Model 1 */

// basic probit
logit absout i.ccp##c.newpow bilat ww2 contig totalmids, cluster(ndyadid) nolog

	
/* Margins Plot for CCP=0 and CCP=1 */

set more off
margins, at(newpow=(0.55(0.05)1)ccp=(0 1) bilat=0 totalmids=1 contig=1)
marginsplot, l(90)


/* Note that the graph is edited within Stata. Color scheme comes from the blindscheme add-on (plottig scheme) */
